Size of Plastic Events in Strained Amorphous Solids at Finite Temperatures 
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We address the system-size dependence of typical plastic flow events when an amorphous solid is 
put under a fixed external strain rate at a finite temperature. For system sizes that are accessible to 
numerical simulations at reasonable strain rates and at low temperatures the magnitude of plastic 
events grows with the system size. We explain however that this must be a finite size effect; for 
larger systems there exist two cross-over length-scales fi and £2, the first determined by the elastic 
time-scale and the second by the thermal energy-scale. For system of linear size L larger than any 
of these scales the magnitude of plastic events must be bounded. For systems of size L > ( there 
must exist (L/£) d uncorrelated plastic events which occur simultaneously. We present a scaling 
theory that culminates with the dependence of the cross-over scales on temperature and strain rate. 
Finally we relate these findings to the temperature and size dependence of the stress fluctuations. 
We comment on the importance of these considerations for theories of elasto-plasticity. 



Introduction: The issues of the statistical correla- 
tions of plastic flow events in strained amorphous solids 
are central to the possible form of the dynamical the- 
ory of elasto-plasticity [1-6]. As such they were at the 
center of extensive research in recent years [7-11]. The 
crucial question is whether these events are spatially lo- 
calized and statistically independent, as assumed often 
in the theoretical development, or are they statistically 
correlated to form extended events that depend on the 
system size. Of particular relevance to the present Letter 
is Ref. [12] in which the authors studied the question for 
zero temperature as a function of the strain rate. At low 
strain rates 7 the plastic events were shown to be spa- 
tially correlated with a system size dependence. At high 
strain rates (compared to elastic relaxation times) the 
correlation were cut-off proportional to 7 _1 / d where d is 
the space dimension. Two crucial questions that remain 
are (i) what is the effect of temperature on this issue. 
Should temperature fluctuations also cut-off the statisti- 
cal correlations? and (ii) if temperature effects do cut off 
the magnitude of plastic flow events, which of the cut-offs 
dominates at a given temperature and strain rate? 

The aim of this Letter is to address these two questions. 
We will show that temperature effects are as important, 
if not more important, in checking the magnitude of plas- 
tic events as the effect of a finite 7. We will present below 
some quantitative estimates of the various effects to com- 
pare their efficacy in bounding the magnitude of plastic 
flow events at a given temperature and strain rate. 

Summary of the Athermal, Quasi-static Simu- 
lations : At athermal conditions T = an amorphous 
solid subjected to very slow strain rate (quasi-static in 
the limit) tends to set up an elasto-plastic steady state 
in which short elastic intervals in which the energy and 
the stress slowly increase are interrupted by plastic flow 
events during which the energy and the stress decrease 
on the short time scale of elastic relaxation. During 
the steady state one can measure accurately the aver- 
age stress drops (Act) or the average energy drops (AU). 
In both two-dimensions [11] and three-dimensions [13] it 



was found that these averages depend on the total num- 
ber of particles as power-laws, 



(AU) = eN a , (Act) = sN 13 



(1) 



with a > and /3 < 0, where e is the mean energy drop 
per particle, and s is a stress scale to be computed below. 
A scaling relation a — (3 = 1 follows from the average 
energy balance equation, cf. [11] 



cty(Act) 



V = (AU) , 



(2) 



where o~y is the flow stress (the mean stress in the ather- 
mal steady state) and fi is the shear modulus. The actual 
values of the exponents a and can depend on the de- 
tails of the inter-particle potential. Typical values of a 
are a bit less than 0.4 in two dimensions [11] and a bit 
more than 0.4 in three dimensions [13]. In Ref. [11] it 
was shown that the number of particles participating in 
a plastic flow events scales like (AU) . 

The effect of finite strain rate: As said in the 
introduction, Ref. [12] showed that finite strain rates 
may cut-off the magnitude of plastic flow events. To 
understand this effect we start by substituting Eq. (1) in 
Eq. (2) to obtain the scale s, 
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Consider next the rate at which work is being done at 
the system and balance it by the energy dissipation in 
the steady state, 



cryjV = (AU)/ Tpl 



(4) 



where r p i is the average time between plastic flow events. 
This time is estimated as the clastic rise time which is 



Tpl 
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(5) 
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We increase our confidence in this estimate by substitut- 
ing it into Eq. (4) together with the other estimates, to 
find perfect consistency. 

Next we note that r p i decreases when N increases. On 
the other hand there exists another crucial time scale in 
the system, which is the elastic relaxation time 



Eq. (9) and then solving for £2 we find 



r cl ~ L/c 



(6) 



where c is the speed of sound c = \J fx/p- Obviously 
this time scale increases with N like N 1 ^. There will 
be therefore a typical scale £1 such that for a system of 
scale L = £1 these times cross. At that size the system 
cannot equilibrate its elastic energy before another event 
is triggered, and multiple avalanches must be occurring 
simultaneously in different parts of the system, each of 
which has a bounded magnitude. We estimate £1 from 
T e i ~ T p i, finding 



(a/ C ) 
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(7) 



Using now the obvious fact that A(£i) ~ (£i/(A)) d we 
compute 
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We first observe the singularity for quasi-static strain 
when 7 — > 0, where £1 tends to infinity, in agreement 
with the results of quasi-static calculations. Thus at low 
temperatures, before the thermal energy scale becomes 
important, the size of plastic flow events can be huge 
indeed. We show next that thermal effects put a much 
more stringent bounds on the magnitude of plastic flow 
events. 

The effect of finite temperatures: The typical 
scale of thermal fluctuations is ksTN where ks is Boltz- 
mann's constant. Comparing with the N dependence of 
the typical energy drop due to plastic flow events, we see 
that the former increases faster with N, and it will catch 
up when 



eN a ~ Nk R T 



(9) 



This equality will hold when the system size L = £2, 
where (^/(A))" 1 = N. Substituting the last equality in 



(A) 



k R T 



(10) 



Recalling that /3 = a — 1 is negative, we again notice 
the singularity at T — > in agreement with the athermal 
quasi-static simulations. 

A Model Glass Example: To put some size esti- 
mates on these crucial length-scales, and to test their 
consequences, we need to choose a model glass. To this 
aim we employ a model system with point particles of 




FIG. 1: A typical equilibrium configuration with 65,536 par- 
ticles. The particles are all point objects, and the ball around 
each particle is of radius \i . 



equal mass m and positions in three-dimensions, in- 
teracting via a pair- wise interaction potentials. In our 
three-dimensional simulations each particle i is assigned 
an interaction parameter A^ from a normal distribution 
with mean (A). The variance is governed by the poly- 
dispersity parameter A = 15% where A 2 = ^ A '^y '- . 

With the definitions = \r.i — Tj\ and Ay = |(A; + Xj), 
the potential assumes the form 



[XiA k fe(fc+2) /BflNltli friA' 
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In our two dimensional simulations below we use the same potential but choose a binary mixture model with 
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a 'large' and a 'small' particles such that All = 1-4, 
Als = 1.18 and A55 = 1.00. Below the units of length, 
energy, mass and temperature are (A), e, m and e/fcs 
where ks is Boltzmann's constant. The time units tq are 
accordingly tq = J (m(A) 2 /e). The motivation of this 
somewhat lengthy form of the potential is to have con- 
tinuous first and second derivatives at the built-in cut- 
1 

off of nj = Ay (k/Bo) k + 2 . In the present simulations 
we chose k = 10, £>o = 0.2. The choice of a quartic 
rather than a quadratic correction term is motivated by 
numerical speed considerations, avoiding the calculation 
of square roots. In the 3D simulations below the mass 
density p = mN/V — 1.3, whereas in 2D p = 0.85. In 
all cases the boundary conditions are periodic. In Fig. 1 
we present a typical 3D equilibrium configuration of the 
system with N = 65536. We measured for this 3D sys- 
tem the shear modulus p = 15.7 and therefore the speed 
of sound is c w 3.5. The value of ay at T = is about 
0.7 and the typical value of CTqo at higher temperatures 
is of the order of 0.5. 

The estimate of £1 depends of course on 7. In our 
3D simulations we have used 7 = 5 x 10~ 5 , and for the 
given values of the speed of sound and of ay we estimate 
£i/(A) ~ 2 x 10 5 which translates to about 10 16 particles. 
Obviously this system size is hugely beyond the capabil- 
ities of molecular simulations. One could in principle 
increase 7, but not beyond a y / (y/ppL) [14]. It therefore 
remains elusive to demonstrate the cross-over due to the 
elastic time-scale in numerical simulations. Nevertheless 
one should remember in any attempt of developing an 
athermal theory of elasto-plasticity that the plastic flow 
events are very large, a fact that cannot be disregarded 
with impunity. 

The cross-over scale due to thermal energies is very 
well within the range of system size available in numerical 
simulations. Making the plausible estimate e e we see 
that already at T = 10~ 3 £2 /(A) is estimated (for j3 = 
—8/15 in 3 dimensions [13]) as £2 ~ 10 2 , which translates 
to just 1 million particles. For T = 10 -2 this estimate 
drops down to about 1000 particles. Thus we expect a 
very rapid cross-over from correlated avalanches to un- 
corrected ones as the temperature rises above 10~ 3 . 

Demonstration of the Thermal Cross-over: A 
very interesting and direct way of demonstrating the 
cross-over due to thermal effects is provided by mea- 
surements of the variance of the stress fluctuations as 
a function of the temperature and the system size. This 
variance is defined by 



(So*) = 



(12) 



where a x is the mean stress in the thermal steady state. 
In Fig. 2 and 3 we display 2D and 3D measurements of 
this quantity which is obtained by averaging the square 
of the microscopic stress fluctuations in long stretches of 
elasto-plastic steady-states of the models described above 
at a fixed 7 = 2.5 x 10~ 5 in 2D and 7 = 5 x 10~ 5 in 3D. 
It is evident that the variance of the stress fluctuations 
decreases as a function of N. Under quasi-static and 
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FIG. 2: The variance of the stress fluctuations as a function 
of the system size N for a 2D system and for various tem- 
peratures. The first power-law (data in squares) is obtained 
under athermal quasi-static conditions where we determine 
for the present model /3 = —0.61, 9 = —0.40. The other plots 
go up in temperature as indicated. The plots are displaced 
by a fixed amount for clarity. Note that the slope decreases 
(becoming more negative) as the temperature increases 



-O-t = 0.001 
-e-T = 0.01 
-q-t = 0.1 
T = 0.3 

T= 0.5 




FIG. 3: The variance of the stress fluctuations as a function 
of the system size N for a 3D system and for various temper- 
atures. The plots are displaced by a fixed amount for clarity. 
Note that the slope decreases (becoming more negative) as 
the temperature increases 



athermal conditions the dependence is a power-law 

(da 2 ) ~ N 2f> , (13) 

where 9 ps —0.4 both in 2D and 3D. One should no- 
tice the difference between the exponent characterizing 
the N dependence of y/ (5a 2 ) and of the athermal mean 
plastic stress drop (Act), in the sense that 9^/3. This 
difference is due to very strong correlations between elas- 
tic increases and plastic drops. At higher temperatures 
the data in Figs. 2 and 3 indicate a clear cross-over to 
independent stress fluctuations in which 



(<5ct 2 



for high temperatures 



(14) 



To capture the temperature and size dependence of the 
variance, and to demonstrate unequivocally the thermal 
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FIG. 4: The scaling function g(x), cf. Eq. (18) for the 3D data 
(upper panel) and the 2D data (lower panel). Note the cross- 
over for x of the order of unity as predicted by Eq. (9). The 
power law decrease at low values of x are in agreement with 
the prediction of £ « 0.33 in both cases. The two black lines 
represent the theoretical prediction for the scaling function 
g(x) for i<1 and for x S> 1. 



cross-over, we first need to separate the thermal from the 
mechanical contributions to (5a 2 ). We write 



(5a 2 ) = (5a 2 ) T + (5a 2 ) 



(15) 



where (5a 2 ) t denotes the thermal contribution which can 
be read from Eq. (10) of Ref. [15], i.e. 



(5a 2 



fiT/V 



(16) 



For the mechanical part we introduce a scaling function 
which exhibits a cross-over according to Eq. (9) . In other 



words, we propose a scaling function g(x) to describe the 
system-size and temperature dependence of the mechan- 
ical part of the variance: 

(5a^)(N,T) = s 2 N 2e g(iNP/k B T) . (17) 
The dimensionless scaling function g(x) must satisfy 



g^; for x - 
gox'' for x 



oo , 
> . 



(18) 



The first of these requirements guarantees that the fluc- 
tuation are in accordance with the athermal limit. The 
second requirements guarantees that after the cross-over 
the fluctuations of the stress become intensive, requiring 
C = -(1 + 26) //3. We compute C ~ 0.33 both in 2D and 
3D. 

We present tests of the scaling function for both our 
2D and 3D simulations in Fig. 4. Examining the scaling 
functions in Figs 4 we see that although the data collapse 
is not perfect, the thermal cross-over is demonstrated 
very well where expected, i.e. at values of a: of the order of 
unity. The asymptotic behavior of the scaling functions 
agrees satisfactorily with the theoretical prediction for 
both the 2D and the 3D data. 

We thus conclude this letter by reiterating that the 
thermal cross-over appears much more aggressive than 
the shear-rate cross-over in cutting off the sub-extensive 
scaling of the shear fluctuations and mean drops. For 
macroscopic systems it should be quite impossible to 
observe plastic events that are correlated over the sys- 
tem size except for extremely low temperatures in the 
nano-Kelvin range. On the other hand nano particles 
of amorphous solids may show at low temperatures and 
low strain rates some rather spectacular correlated plas- 
tic events. 
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